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Abstract 

We study the spatial correlator of the topological charge density 
operator in pure SU(3) gauge theory and in two flavor QCD. We show 
that the data for distances up to about 1 fm is consistent with a vacuum 
consisting of individual instantons and closely bound pairs. The 
percentage of paired objects is twice as large on the dynamical 
configurations than on the pure gauge ones, implying increased molecule 
formations due to fermionic interactions. 



Our understanding of the QCD vacuum has increased considerably in the 
last few years. Lattice calculations using different algorithms have consistent 
predictions for the topological susceptibility in pure gauge QCD G|, and the 
first results in dynamical systems have been published recently 0,^- Results 
concerning the density and size distribution of instantons are less consistent, 
but there is increasing evidence that the pure gauge QCD vacuum is filled with 
instantons of average radius ~ 0.3 fm, with a density of about 1 fm~ 4 

Instanton Liquid Models (ILM) describe the phenomenological properties of 
instantons in the QCD vacuum closely Q. Even the Random ILM provides 
an accurate description of the pure gauge instanton vacuum indicating that the 
gauge interaction and consequently the spatial correlation between instantons is 
small. The situation is quite different for systems with dynamical light quarks. 
In the zero quark mass limit the topological susceptibility is zero and there are 
no unpaired topological objects in the vacuum. Yet one expects instantons to 
be present. Chiral symmetry is spontaneously broken, < iptp >^ in the zero 
quark mass limit. According to the Casher-Banks formula, the chiral condensate 
(summed over quark flavors) is 



< >= 7rp(0), (1) 

implying that yo(0), the density of the eigenmodes of the Dirac operator at 
eigenvalue zero, is finite. Since instantons are the leading candidates to create 
near-zero eigenmodes of the Dirac operator, the vacuum is likely to be filled 
with instantons. Fermions create an attractive force between oppositely charged 
instantons. It is therefore natural to expect that the instantons of the zero quark 
mass vacuum form instanton-antiinstanton pairs, molecules. In the case of finite 
quark mass the topological susceptibility does not vanish but it is proportional 
to the quark mass || 

X =< W > ^ + 0(mi) = 0(mi), (2) 

where f v = 132MeV is the pion decay constant. The vacuum in this case has 
unpaired objects in addition to the molecules. 

Even though the above picture describing spatial correlation between instan- 
tons sounds very natural, no evidence from lattice calculations has supported it 
so far. The culprit is most likely the lattice approach. In order to reveal individ- 
ual topological objects, vacuum fluctuations have to be removed, the vacuum 
has to be smoothed. Almost all smoothing procedures distort the vacuum, they 
destroy molecules especially easily. In addition the smoothed configurations are 
frequently analyzed by a pattern-recognition algorithm to identify individual 
instantons. Most algorithms have a built-in cut-off which limits the nearest 
objects it can resolve, further limiting the possibility of finding closely bound 
pairs. 
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In this paper we study the topological density correlator 

C(r) = i / d 4 xq(x)q(x + z), r = \z\, (3) 

where q(x) is the topological charge density measured using the improved charge 
operator of Ref . || , |9) . We measured C (r) on two sets of configurations with 
similar lattice spacings. The first set is a pure gauge ensemble of 16 3 x 32 con- 
figurations at /3 = 6.0, the is a dynamical ensemble of the same size, generated 
with two flavors of staggered fermions at j3 = 5.7 with ma = 0.01. The lattice 
spacing of the pure gauge ensemble is a = 0.095 fm, of the dynamical ensemble 
is a = 0.11 fm. Both ensembles are from the NERSC QCD archive [jro|. For- 
tius study 170 of the pure gauge configurations [jll| and 83 of the dynamical 
configurations jl2| (the entire available data set) were used. 

The gauge configurations have to be smoothed to reveal the topological 
content. Here we used APE smearing |u| with parameter c = 0.45 as the 
properties of this smoothing has been extensively tested in Refs. ra^Pt- N 
steps of APE smearing smoothes a configuration to distance d s ~ Q 
but preserves the properties of the vacuum at longer distances. First we will 
compare the topological density correlators on the pure gauge and dynamical 
configurations after the same number, N — 30, APE steps. Since the lattice 
spacings of the two ensembles are approximately the same, that corresponds 
to about the same physical smoothing and even observables that change with 
smoothing can be compared. Afterwards we will vary N between 10 and 60 to 
justify the above assumption. 

First we consider the topological susceptibility. The susceptibility is largely 
independent of the number of smoothing steps on both configuration sets. Sev- 
eral previous studies found that the topological charge has very large autocor- 
relation time on dynamical configurations, sometimes hundreds of molecular 
dynamics time units. We did not find this problem here. The average charge 
on the dynamical configurations is < Q >= 0.16 ± 0.25 and the charge distri- 
bution also appears standard. We found x 1 ^ 4 = 193(4)MeV on the pure gauge 
ensemble and x 1 / 4 = 130(5)MeV on the dynamical ensemble. This decrease 
of the susceptibility is expected according to Eq. |^. Using the published pion 
mass value m^a = 0.25 of the dynamical ensemble Jl2| in Eq. we obtain 
f n = (105 ± 20) MeV, fairly close to the experimental value f v = 132 MeV. 
(Note that Eq. || differs from the equivalent equation of Ref. Q []) 

Figures [j] and ^ show C(r) versus r for both configuration sets after N = 30 
APE smoothing steps. The inserts of the figures enlarge the large r tail. At this 
smoothing level the vacuum fluctuations are largely removed. A few properties 
of the instanton vacuum can be immediately deduced from the figures. At small 
distances C (r) is dominated by the auto-correlator of individual instantons and 

lr The author is indebted to Peter Hasenfratz for checking Eq. El 
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Figure 1: Spatial correlation of the topological charge density on pure gauge 
configurations. The insert enlarges the large r region. The solid line corresponds 
to the four parameter fit described in the text. 
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Figure 2: Same as Figure |l| but for the dynamical ensemble. 



Table 1: The results of the four parameter fit for the pure gauge and the dy- 
namical configurations. 





Pure gauge 


Dynamical 


N /V [fur 4 ] 


1.08±0.05 


0.35± 0.02 


p[fm] 


0.25±0.01 


0.30± 0.01 


d[fm] 


0.53± 0.03 


0.61± 0.02 


2n p /N 


0.16± 0.02 


0.29± 0.03 



C(r) can be approximated by a dilute gas. In the dilute gas approximation 
the height of the correlator is proportional to the number of topological objects 
of the configuration, iVo, and the width is related to the average radius of the 
instantons, p. At distances r ~ 2p, the correlator probes the close neighbors 
of the instantons. If there is no spatial correlation between topological objects, 
C(r) should be about zero for r > 2pf\. If pairing of oppositely charged objects 
occur, C(r) is expected to be negative, reaching its minimum around r ~ 2p. 
Comparing now the pure gauge and dynamical correlators of figures |l| and |^ we 
observe that the widths of the correlators are almost identical but the height at 
r = is very different. This implies that p is approximately the same for the 
two ensembles but the dynamical configurations have about three times fewer 
topological objects. The tails of the two distributions reveal further differences. 
While C(r) on the pure gauge ensemble is consistent with zero for r > 0.7 fm, 
C(r) on the dynamical ensemble shows an unmistakable negative dip signaling 
opposite-charge correlation. 

To quantify these differences, we model the vacuum with a simple picture. 
Assume that there are Nq topological objects in the vacuum with uniform radius 
p. If the topological density of a single instanton of radius p centered at the origin 
is qo(x) and the instantons are non-overlapping, the topological density of the 
configuration in this model is 

N 

= X! s i1o( x - Xi), (4) 

where Xi is the center of the the ith topological object and s, = ±1 depending 
on whether it is an instanton or antiinstanton. The topological correlator is 

C(z) = SjSjC p (z - Xj+ xj) (5) 

where C p (z) is the topological correlator of a single instanton. If the instantons 
are randomly distributed, averaging over configurations will give 

<C{z) >= NoC p (z). (6) 

2 In the continuum reflexion positivity requires C(r) < for all r ^ 0. On the lattice that 
is not the case for small r. Since smoothing removes most of the vacuum fluctuations one 
expects C(r) to be dominated by non-perturbative vacuum structures. 
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If only N — n p objects are randomly distributed and the other n p objects form 
molecules (i.e. 2n p oppositely charged instantons are paired), the average of the 
topological correlator is 

< C{z) >= N C p {z) - n p C d (z), (7) 

where Cd(z) is the correlator of an instanton-instanton pair separated by dis- 
tance d, averaged over the direction of d. C p and Cd can be calculated using 
the analytical form of the single instanton charge distribution. If the instanton 
size is p ~ 0.3 fm and the pairs are closely bound, i.e. d ~ 0.6 fm, this model 
can be valid in the < r < 0.6 — 0.9 fm range. The measured C(r) then can 
be fitted with the four parameters Nq, n pi p and d. In the fit we keep data 
with 0<r/a<5 — 8asforr/a>L/2 = 8 finite size effects become important. 
Varying the upper range of the fit between r/a — 5 and 8 hardly effects the 
fit parameters, it changes the x 2 OI the fit only. The solid lines in figures [l] 
and H correspond to the best fit. The fit parameters listed in table |l| confirm 
our earlier observations. The average instanton size on both ensembles is about 
0.3 fm. The instanton density on the pure gauge ensemble is about 1 fm -4 , in 
agreement with expectations, while the density on the dynamical configurations 
is about a third of that. On both ensembles we observe pairs at a distance of 
about 2p. On the pure gauge configurations 15 percent of the instantons are 
in pairs, while on the dynamical configurations close to 30 percent are paired 
up. It is not surprising to find pairs even on the pure gauge configurations, as 
there is an attractive gauge interaction between oppositely charged instantons, 
but pairing is clearly enhanced on the dynamical ensemble. 

Until now we have compared the dynamical and pure gauge ensembles at 
the same level of APE smoothing. In the remainder of this paper we study the 
dependence on the level of smoothing. The smoothing has two different effects 
on a configuration: it removes vacuum fluctuations and it also distorts and 
annihilates instantons. These two effects cannot be separated. The first effect 
is the important one for small amounts of smoothing, but as the configuration 
is smoothed many times, the second one becomes more relevant. We have 
measured the correlator C(r) at N = 10,20,30,40,50 and 60 levels of APE 
blocking on both ensembles and also at N = 15 on the dynamical configurations 
only. The same level of APE smoothing can have very different effect at different 
lattice spacing. The smearing distance d s = ay/Nc/3 characterizes the removal 
of vacuum fluctuations but it is not clear what is the variable (if one exists at 
all) that describes how instantons are annihilated or otherwise distorted by APE 
smoothing at different lattice spacings. In the following we report the results as 
the function of the smearing distance. The qualitative form of the correlator is 
the same at all levels of smearing and we fit the data with four parameters as 
described above. The fit is excellent in every case except at the N = 60 level 
smoothing where the smearing distance is comparable to the lattice size and 
finite size effects become important. 

Figure shows the ratio of paired objects and Ao as the function of d s . 
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Figure 3: 2n p /No as the function of smearing distance for the dynamical (dia- 
monds) and pure gauge (crosses) ensembles. 



The ratio decreases rapidly up to d s ~ 0.2 fm but changes slower after that, 
indicating that most of the vacuum fluctuations have been removed and further 
change is due to the annihilation of topological objects. Comparing data at 
the same smoothing distance even increases the previously observed difference 
between the pure gauge and dynamical ensembles; the dynamical configurations 
show more than twice the pairing rate of the pure gauge ones. 

It is interesting to ask, which instantons contribute to the topological sus- 
ceptibility, all of them or only the unpaired ones? Figure || shows the density of 
the unpaired objects, (No — 2n p )/V as a function of the smearing distance d s . 
The horizontal lines indicate the value of the corresponding topological suscep- 
tibilities in fm . The agreement implies that the unpaired instantons follow a 
Poisson distribution and only they contribute to the topological susceptibility. 

We have studied only one set of dynamical ensemble with two quark flavors. 
It would be important to extend this work to include different, preferably smaller 
quark masses and four quark flavors. In both cases one expects the fermionic 
effects to be stronger. We mention here that the NERSC archive contains 
two smaller (less then 50 configurations) dynamical configuration sets, both 
at heavier quark masses. We have analyzed the available configurations, but 
within errors they show no significant difference compared to the dynamical 



7 



1.5 



- £ 



1.0 



X X 



I 



0.5 



0.0 



0.1 



0.2 0.3 
djfm] 



0.4 



Figure 4: (No — 2n p )/V as the function of the smearing distance for the dy- 
namical (diamonds) and pure gauge (crosses) ensembles. The horizontal lines 
indicate the value of the corresponding topological susceptibilities. 



configurations discussed here. 

In summary, we have demonstrated important differences between pure 
gauge and dynamical QCD configurations. On the dynamical configurations 
the topological susceptibility is considerably smaller than on the pure gauge 
configurations. This is due both to the decreased density of instantons and to 
the stronger pairing of oppositely charged objects. In the case we studied here 
we found that the density of instantons decreased by about a factor of three 
while pairing occurred twice as frequently on the dynamical configurations. 

The author is indebted to T. DeGrand for carefully reading the manuscript. 
This calculation would not have been possible without the configurations avail- 
able at NERSC. We would like to thank the Colorado High Energy experimental 
group for allowing us to use their work stations. This work was supported by 
the U.S. Department of Energy. 
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